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Abstract 

An intensive spectroscopic study was performed for three representative solar twins (HIP 56948, 
HIP 79672, and HIP 100963) as well as for the Sun (Moon; reference standard), with an intention of 
(1) quantitatively discussing the relative-to-Sun similarities based on the precisely established differential 
parameters and (2) investigating the reason causing the Li abundance differences despite their similarities. 
It was concluded that HIP 56948 most resembles the Sun in every respect including the Li abundance 
(though not perfectly similar) among the three and deserves the name of "closest-ever solar twin" , while 
HIP 79672 and HIP 100963 have somewhat higher effective temperature and appreciably higher surface Li 
composition. While there is an indication of Li being rotation-dependent because the projected rotation in 
HIP 56948 (and the Sun) is slightly lower than the other two, the rotational difference alone does not seem 
to be so large as to efficiently produce the marked change in Li. Rather, this may be more likely to be 
attributed (at least partly) to the slight difference in T e s via some T e ff-scnsitivc Li-controlling mechanism. 
Since the abundance of beryllium was found to be essentially solar for all stars irrespective of Li, any 
physical process causing the Li diversity should work only on Li without affecting Be. 

Key words: stars: abundances — stars: atmospheres — stars: individual (HIP 56948, HIP 79672, 
HIP 100963) — stars: rotation — stars: solar analog 



1. Introduction 

Can we find such a star that indiscernibly resembles 
our Sun in every respect? This "solar twin 1 survey" , an 
ever-attracting subject for stellar astronomers, has made 
significant progress since 1990s, thanks to the improve- 
ment in the precision of stellar parameter determinations. 

Since Porto de Mello and da Silva (1997) reported the 
remarkable similarity of HIP 79672 (= 18 Sco = HR 6060 
= HD 146233; V = 5.50) to the Sun, this star has main- 
tained the status of best solar twin candidate almost for a 
decade (see also Soubiran & Triaud 2004). In the mean- 
time, by using a numerical technique developed by Takeda 
(2005; hereinafter referred to as Paper I) for establishing 
the parameter differences between two similar stars with 
high precision, Takeda et al. (2007; hereinafter Paper 

* Based on data collected at Subaru Telescope, which is 
operated by the National Astronomical Observatory of 
Japan. The electronic tables El and E2 will be avail- 
able at the E-PASJ web site upon publication, while 
they are provisionally placed at the WWW site of 
(http:/ /optik2. mtk.nao.ac.jp/ takeda/solartwins/). 

1 We use the term "solar twin" for those special solar-type stars 
which have particularly high similarity to the Sun with respect 
to spectra as well as stellar parameters. See Appendix A of 
Takeda et al. (2007) and the references therein for the literature 
concerning this theme. 



II) conducted a comprehensive study of solar analog stars 
and found that HIP 100963 (= HD 195934; V = 7.09) is 
an equally good (or even better) solar twin as HIP 79672. 

Yet, there is one concern. While these HIP 79672 and 
HIP 100963 are certainly very similar to the Sun in terms 
of the stellar parameters and the general appearance of 
the spectra, one marked dissimilarity exists in a particular 
part of the spectrum: the strength of Li line at 6707.8 A in 
these two stars is appreciably stronger compared to the so- 
lar case (cf. figure 3 in Soubiran & Triaud 2004 and figure 
A. 2 in Paper II). This decisive difference in the surface 
Li abundance is actually "a fly in the ointment," which 
makes us somewhat hesitate to regard them as "real" so- 
lar twins. 

Interestingly, however, Mclendcz and Ramirez (2007) 
recently reported that HIP 56948 (= HD 101364; V = 
8.70) appears to be an ideal solar twin in the sense that it 
has essentially solar parameters and the low Li abundance 
similarly to the Sun. If this is confirmed, this star may 
deserve being called as a genuine solar twin. It would 
thus be worth carrying out an independent check analysis 
in order to ascertain whether HIP 56948 really resembles 
our Sun on every point including the Li abundance. 

Another related subject of interest is the cause of such a 
difference in the Li abundance among these superficially 
very similar solar twins. It was concluded in Paper II 
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based on the analysis of 118 solar- analog dwarfs around 
carly-G type that the surface Li abundance is closely cor- 
related with the macroscopic line-broadening parameter 
(comprising macroturbulcnce plus rotation); i.e., the sur- 
face Li tends to be higher (i.e, less depleted) as the line- 
width becomes broader (cf. figure 13 therein). Since the 
macroturbulence (due to the granular motion of stellar 
convection origin) is unlikely to differ much among simi- 
lar solar-type stars, this observational fact suggests that 
the most decisive factor controlling the surface Li abun- 
dance is the stellar rotation or the angular momentum 
(i.e., faster rotation tends to suppress the envelope mix- 
ing and leads to less depletion of Li). Then, is the distinc- 
tion between the Li-strong (HIP 79672 and HIP 100963) 
and Li-weak (HIP 56948 and naturally the Sun itself) so- 
lar twins simply caused by the difference in the rotational 
velocity? This point should be checked by careful deter- 
minations of the projected rotational velocities of these 
stars. 

Besides, we should also pay attention to two other 
related viewpoints in connection with this "rotation- 
mixing-surface Li" relationship. The first is the stellar 
activity which tends to diminish/enhance as the rotation 
becomes slower/faster. If the rotation is the key factor af- 
fecting the surface Li, does the Li-strong solar twins show 
higher activity than Li-weak ones? The second is the sur- 
face abundance of beryllium, which is destroyed when con- 
veyed into the hot stellar interior by envelope mixing sim- 
ilarly to lithium at temperature of T ~ 3.5 x 10 6 K (higher 
than the case of Li which is burned at T ~ 2.5 x 10 6 K). 
It is interesting to see whether any difference is observed 
in the surface abundance of Be between the Li-strong and 
Li- weak groups, which may provide us with an observa- 
tional constraint on the origin of Li discrepancies among 
these solar twins 

Motivated by these considerations, we decided to con- 
duct an intensive spectroscopic study for these representa- 
tive solar twins (HIP 56948, HIP 79672, and HIP 100963 
along with the Sun/Moon as the comparison standard) 
based on the high-dispersion spectra obtained by the 
Subaru Telescope with HDS, in order to (1) quantitatively 
discuss the relative-to-Sun similarities of these three stars 
while precisely establishing their "star— Sun" differential 
parameters by applying the method of Paper I and to (2) 
investigate the reason/mechanism causing the difference 
between the Li-strong and Li-weak groups by examining 
the rotational velocity, the degree of stellar activity, and 
the Li as well as Be abundance. This is the purpose of 
this study. 

2. Observational Data 

The observations of HIP 56948, HIP 79672, HIP 100963 
and the Moon (substitute for the Sun) were carried out 
in the night of 2008 June 15 (Hawaii Standard Time) by 
using the High Dispersion Spectrograph (HDS; Noguchi 
et al. 2002) placed at the Nasmyth platform of the 8.2-m 
Subaru Telescope, which can record high-dispersion spec- 
tra covering a wavelength portion of ~ 1600 A (blue cross 



disperser) or ~ 2600A (red cross disperser) with two CCDs 
of 2Kx4K pixels at a time. 

In order to cover the wide wavelength range from near- 
UV (~ 3000A) to red (~ 7000A), each star was ob- 
served at two different wavelength settings (standard Ub 
with blue cross disperser for ~ 3000-4500 A and stan- 
dard Yc with red cross disperser for ~ 4400-7000 A). 
With the slit width set at 0."4 (200 fjxa) and no on- 
chip binning of pixels, the resolving power of the ob- 
tained spectra is R ~ 90000. The integrated exposure 
times for each star at Ub/Yc settings are 64 min/32 min, 
18 min/7 min, 32 min/ 16 min, and 2 min/0.5 min for 
HIP 56948, HIP 79672, HIP 100963 and the Moon, re- 
spectively. 

The reduction of the spectra (bias subtraction, flat- 
fielding, scattered-light subtraction, spectrum extraction, 
wavelength calibration, continuum normalization) was 
performed by using the "cchellc" package of the software 
IRAF 2 in a standard manner. The estimated S/N ra- 
tios at each of the wavelengths calculated as the square 
root of the resulting photoelectron counts (ADN x gain) 
are graphically depicted in figure 1. We can sec from 
this figure that sufficiently high S/N ratios of ~ 500-1000 
are achieved in the most sensitive red region, though this 
value is considerably reduced even by a factor of ~ 10 at 
- 3100 A of near-UV where Be n lines arc located. 

3. Parameter Determination 

3.1. Standard Stellar Parameters 

Following the procedure described in subsubsection 
3.1.1 of Paper II, we measured the equivalent widths 
(EW) of Fe I and Fe II lines on the "Yc setting" spec- 
tra covering ~ 4400-7000 A. Based on these EW values, 
the four standard 3 atmospheric parameters [T^ d (effec- 
tive temperature), \ogg (surface gravity), vt (microtur- 
bulent velocity dispersion), and {Fe/H} std [= Af c d - 7.50] 
(Fe abundance )], which are necessary for constructing 
model atmospheres, were spectroscopically derived by us- 
ing the TGVIT program (Takeda ct al. 2005; cf. section 
2 therein). This method is based on the principle search- 
ing for the most optimum solution in the 3-dimcnsional 
(T c ff, log g, v t ) space such that simultaneously satisfy- 
ing the three requirements of (i) the excitation equilib- 
rium, (ii) the ionization equilibrium, and (iii) the EW- 
independence of the abundances (cf. Takeda, Ohkubo, & 
Sadakane 2002). The detailed EW data and the Fe abun- 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 

3 We use the notation of "standard" here, which means that these 
have characters of "absolute" parameters in the usual sense, in 
order to discriminate them from the "differential" parameters 
(relative to the Sun) presented in subsection 3.5. 

4 We intentionally expressed this quantity as {Fe/H}, in order 
to clarify that it is still an absolute quantity (i.e., essentially 
equivalent to ^4p^ d ) and should be distinguished from the strictly 
differential metallicity relative to the Sun, [Fe/H] (= Ap e ar — 
A F " n ), which is also denoted as AAp c in subsection 3.5. 
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Fig. 2. Comparison of the equivalent widths of Fc I and Fe II 
lines measured on the Subaru/HDS spectra for the determi- 
nations of atmospheric parameters in this study (abscissa) 
with those used in the previous investigations based on the 
OAO/HIDES spectra (ordinate), (a) HIP 79672 (Paper II), 
(b) HIP 100963 (Paper II), (c) Moon (Paper II), and (d) Moon 
(Paper I). 



Fig. 1. Distribution of S/N ratios (estimated as the square 
root of photoelcctron counts) of the spectra used for this 
study, which are divided into four wavelength regions: 
~ 3000-3700 A (blue CCD) and ~ 3700-4600 A (red 
CCD) in the Ub setting; ~ 4400-5700 A (blue CCD) and 
~ 5700-7000 A (red CCD) in the Yc setting. The absence 
of data at ~ 3700A and ~ 5700A corresponds to the joint 
of mosaicked CCDs. Spurious spikes seen at several wave- 
lengths are due to bad columns of CCDs, (a) HIP 56948, (b) 
HIP 79672, (c) HIP 100963, and (d) Moon. 

dances corresponding to the final parameters for each star 
are presented in electronic table El. The results are sum- 
marized in table 1, where the related stellar parameters 
(L, M, and age) evaluated as in subsection 3.3 of Paper 
II are also given. 

Comparing the present EW data measured from the 
Subaru/HDS spectra (i? ~ 90000, S/N - 500-1000) with 
those from OAO/HIDES spectra in Paper II and Paper 
I {R ~ 70000, S/N - 200-600), we confirmed a general 
consistency as shown in figure 2. However, a close inspec- 
tion revealed a slight systematic difference in the sense 
that EW(Subaru) tends to be by ~ 1-2% smaller than 
EW(OAO). According to this delicate systematic change, 
marginal differences are seen in the present results of such 
absolute parameters when compared to those in these pre- 
vious papers; e.g., for the case of the Sun/Moon, T e g and 
{Fc/H} have been lowered by ~ 30 K and ~ 0.03 dex, re- 
spectively, though these changes should not matter in the 
differential analysis (subsection 3.5). 

3.2. Rotational Velocity 

In order to evaluate the stellar projected rotational 
velocity (u e sini), we determined the total macrobroad- 
cning parameter, vm, which is the e-folding width 



of the Gaussian macrobroadening function, /m(^) oc 
exp[— (v/vm) 2 }, by way of the line-profile fitting as done 
in subsection 4.2 of Paper II. Unlike the previous case 
(where the fitting was applied to the spectrum portion 
at the 6080-6089 A region), however, we performed the 
fitting analysis to each of the "individual" Fe I and Fe n 
lines (the same lines as used for the EW measurements 
in subsection 3.1) as Takeda (1995) did for the solar flux 
spectrum, because the macroturbulence (to be subtracted 
from the total macrobroadening) is considered to be dif- 
ferent from line to line because of its depth-dependence 
(cf. Takeda 1995). 

We used the line-broadening model adopted by Takeda 
ct al. (2008). That is, the total macrobroadening func- 
tion, /m(v), is assumed to be the convolution of three 
Gaussian component functions f a oc exp[— (v/v a ) 2 }, where 
a is any of "ip" (instrumental profile), "rt" (rotation), and 
"mt" (macroturbulence); i.e., 

2 2,2.2 / 2,2 \ /,i 

V M = v ip + «rt + V (= v ip + u r+m) . I 1 ) 

where v r+m is the "macroturbulcnce+rotation" parame- 
ter used in Paper II. These broadening parameters (v- lp , 
v T t, and v m t) may be related to the more realistic quan- 
tities as Vi p ~ (c/i?)/(2\/ln2) (2.00 km s _1 in the present 
case of R ~ 90000), i> r t — 0.94u e siru (v e and i are the 
equatorial rotation velocity and the inclination angle), 
and v m t ~ 0.42£rt (Crt: radial-tangential macroturbu- 
lence dispersion; cf. Gray 2005), as explained in footnotes 
10 and 12 of Takeda et al. (2008). 

Further, since we may reasonably postulate that the 
macroturbulence velocity field in the solar atmosphere can 
be applied to all of the three solar twin targets, we assume 
an analytical form of the depth-dependent macroturbu- 
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Fig. 3. Line-broadening parameters (derived by the pro- 
file-fitting method applied to blend-free Fe lines) plotted 
against the mean-depth of line formation ( (log T5000 ) ) 1 where 
green open circles, blue filled circles, and red dots rep- 
resent (rotation+macroturbulence), v m t (macroturbu- 
lent velocity parameter, supposed to be depth-dependent as 
Vmt = 1-60 — 0.11 log rsooo — 0.19 log r| 000 ; cf. subsection 
3.2) and v T t (rotational broadening parameter obtained as 
i v i+ m ""mt) 1 ' 2 )' respectively. The blue horizontal line indi- 
cates the average value of v v t {{v T t); cf. table 1) calculated for 
the lines of (lognsooo) < -0.7. (a) HIP 56948, (b) HIP 79672, 
(c) HIP 100963, and (d) Moon. 

lence (in terms of T5000, the optical depth at 5000 A) 

w m t(=0.42C RT ) = 1.60-0.111ogT 50 oo-0.191ogT 5 2 000 ,(2) 

since £rt may be approximately expressed as 3.8 — 
0.251o g7"5ooo~ 0.451ogr| 000 from figure 2 of Takeda (1995). 
Now the procedure makes as follows. 

(1) First, vm is determined by applying the profile-fitting 
method to a blend free Fe line. 

(2) Then, i> r + m (macroturbulence plus rotation) is ob- 
tained by subtracting the instrumental broadening effect 

(vi p = 2.00 km s _1 ) as v t+m = yjv^-vfp. 

(3) From the measured equivalent width, the mean-depth 
of line formation ((logr 50 oo)) relevant for this line can be 
computed (cf. subsection 5.1 in Takeda 1995), which is 
sufficient to assign an appropriate value of v m t to this line 
with the help of equation (2). 

(4) Eventually, v T t is evaluated as v T t = J i'r+ m — w„ t • 
While the detailed results of v r+m and v Tt (along with 

the assigned (logrsooo)) for each line are presented in 
electronic table E2, these w r + m , v m t, and v r t are plotted 
against T5000 m figure 3. We can see from this figure that 
the depth-dependent tendency of t> r + m is almost removed 
in v T t by subtracting the effect of v m t. Finally, we obtained 
(v T t) as the parameter representing w sini 5 by averaging 
the u rt 's with (logrsooo) < —0.7 (deep-forming lines with 



Admittedly, we can not hope to relate the "exact" value of v c s'mi 
to (v r t ) within the framework of such a rough modeling of line- 
broadening functions (all assumed to be the Gaussian form). 
However, we may reasonably expect that (v^t) is proportional 
to v e sin i with a factor not much different from unity. This 
is sufficient for our present purpose, because what we want to 
know is the "differential" characteristics (i.e., the ratio of (i>rt) 
between two stars is considered to be the ratio of actual u e sini). 
At any rate, it is encouraging that the resulting (v r t) value of 
2.13 km s -1 for the Sun/Moon is quite close to the actual solar 
t! e sini value of 1.9 km s -1 , by which we may regard that our 
approximation (suggesting v T t ~ 0.94u e sini) is not bad. 



(logr 5 ooo) > —0.7 were not used for the averaging because 
of the larger uncertainties due to the weakness of the line- 
strength), as given in table 1. 

3.3. Li Abundance 

The portion of the observed spectrum (6706.3- 
6709.3 A) comprising the Li 1 resonance doublet at 
~ 6707.8 A, along with Kurucz ct al.'s (1984) solar flux 
spectrum, is shown in figure 4 (left panels). We can rec- 
ognize from this figure that the strengths of the Li line 
for HIP 79672 and HIP 100963 are markedly larger than 
those for HIP 56948 and the Sun/Moon, classifying these 
four into Li-strong and Li- weak groups. As in Paper II, 
the Li abundance (A-^i) was determined from the Li 1 
doublet lines at ~ 6707. 8A in the same manner as de- 
scribed in Takeda and Kawanomoto (2005). Namely, we 
first establish the LTE abundance (A^™) by applying the 
method of synthetic profile fitting to the spectrum fea- 
ture of Fe 1 + Li 1 lines (see the right panels in figure 
4). Then, the EW(Li 1 6708) is inversely calculated from 
such obtained -A^j™- Finally, while taking into account 
the non-LTE effect, A^ TE is calculated from EW(Li 1 
6708). The resulting for each star is presented in 

table 1. In all the four cases studied, the non-LTE cor- 
rection A(= A^lte _ a lt e ) turned out to be +0.07. The 
solar Li abundance of 0.91 derived in this study based 
on the spectrum of Moon (Subaru/HDS) is in excellent 
agreement with the result of 0.92 concluded by Takeda 
and Kawanomoto (2005) based on the spectrum of Moon 
(OAO/HIDES) as well as the solar flux spectrum (Kurucz 
et al. 1984). 

3.4- Be Abundance 

The spectrum portion of 3129.5-3131.5 A comprising 
Be 11 lines at 3130.42 A and 3131.07 A is shown in figure 
5, where each stellar spectrum is compared with Kurucz 
et al.'s (1984) solar flux spectrum. A glance of this figure 
suffices us to convince that Be line features are essentially 
the same as the solar case for all stars. According to the 
theoretical simulation shown in the lowest panel of this 
figure, the agreement of Ab c with the solar value appears 
to be very good, presumably to within ~ 0.1 dex (though 
the uncertainty may be somewhat larger for HIP 56948 
where the spectrum quality is comparatively poor) . This 
fact clearly suggests that Be makes a clear distinction from 
Li (showing an appreciable difference from star to star) 
in spite of their rather similar characters comparatively 
easily destroyed in the stellar interior, at least for these 
solar twin stars are concerned. This result is consistent 
with what Randich et al. (2002) concluded for carly-G 
dwarf stars in open clusters. 

3.5. Differential Analysis 

Now that the "standard" atmospheric parameters were 
established in subsection 3.1, we can derive the "differ- 
ential" parameters Api^j (p is any of T c g, logg, v t , and 
A-pc) of star i relative to any other arbitrary comparison 
star j by using the method described in Paper I, where 
several practical quantities were defined such as (i) the 
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Fig. 4. Left side: Observed spectra of target stars in the 
spectrum region comprising Li I 6707.8 line (open symbols) 
in comparison with the solar flux spectrum (red dashed 
line). Right side: Theoretically simulated spectrum in the 
Li I 6707.8 line region (blue solid line) fitted with the ob- 
served stellar spectrum (open symbols), (a) HIP 56948, (b) 
HIP 79672, (c) HIP 100963, and (d) Moon. 

average of the direct solution, (Apij) (average of Api-j 
and —Apj—i), (ii) the intermediary solution via star k, 
(ApuQj) (= (Ap ik ) + (Apik)), and (iii) the average of 
the intermediary solution, ((ApiQj)) (average of (Ap^^j) 
over various k). 

Since we are interested in the parameter differences rel- 
ative to the Sun, we take i = 1,2,3 and j = (see table 1 for 
the numbering of each star), and two intermediary stars 
can be assigned for any pair (e.g., for the case of i = 1 and 
j = 0, we can take k = 2 or k = 3). The detailed results 
for HIP 56948 (i = 1), HIP 79672 (i = 2), HIP 100963 
(i = 3) are presented in tables 2, 3, and 4, respectively. 
(Note that these three tables arc formatted in the same 
manner as in tables 2-9 of Paper I.) The values of (Apio) 
(direct solution) and ((ApjQ )) (average of the intermedi- 
ary solution) are separately summarized in table 5, where 
the differences in w rt and Au are also given. It is worth 
noting that the comparison of these two direct and inter- 
mediary solutions may provide us with an opportunity of 
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Fig. 5. Observed spectra of target stars in the near UV re- 
gion comprising two Be II lines at 3130.42 and 3130.07 A (sym- 
bols: filled and open circles correspond the spectra of different 
echelle orders 189 and 190, respectively), in comparison with 
Kurucz et al.'s (1984) solar flux spectrum (red dashed line). 
The continuum level of each spectrum has been so adjusted as 
to achieve a consistency between the stellar and the reference 
solar spectrum. (a) HIP 56948, (b) HIP 79672, (c) HIP 100963, 
and (d) Moon. In the lowest panel (e), theoretically synthe- 
sized spectra for Sun (blue solid lines), are compared with the 
solar flux spectrum (red dashed line), which were computed 
by using the atomic data given in Primas et al. (1997) with 
different Be abundances of [Be/H] = -0.6, -0.3, 0.0, +0.3, 
and +0.6. 
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checking/estimating the accuracy of the results. 
4. Discussion 

4-1. Which Is the Best Solar Twin? 

When we compare the parameter differences (relative 
to the Sun) for HIP 79672 and HIP 100963 given in ta- 
ble 5 with those of Paper II (see table A.l therein), we 
see a notable discrepancy in the values of AT s, despite 
that other Alogg, Avt, and AAp e are mostly in good 
agreement. Namely, while Paper II derived very small 
AT eff (+1.7 K/-1.6 K for HIP 79672/100963), the present 
results (+48.5 K/+38.2 K for the direct solution) indi- 
cate appreciably larger values by ~ 40-50 K than these. 
Although the reason for this AT c ff discrepancy is not clear, 
it would presumably be ascribed to the difference in the 
used EW data set. At any rate, because of the reasonable 
consistency between the direct and intermediary solutions 
(compare the values in the upper and lower rows in tabic 
5) and the use of the spectrum data of much higher quality 
(in terms of both the S/N ratio and the spectrum resolv- 
ing power), we would place larger weight in the present 
results. 

By inspecting table 5, we can summarize as follows con- 
cerning the similarity or dissimilarity of these three pro- 
gram stars to the Sun in terms of each checkpoint. (In the 
discussion of the parameter differences given below, we re- 
fer to the averaged values of the direct and intermediary 
solutions for convenience.) 

• ATcff: HIP 56948 is manifestly more solar like (~ 
+10+10 K) than other HIP 79672 and HIP 100963 
(~ +40 + 10 K). 

• Alogg: All three stars do not show appreciable dif- 
ferences from the solar value, if we consider the nom- 
inal uncertainty of ~ 0.01 dex. 

• Avt - - HIP 79672 shows slightly higher v t than 
the Sun by +0.03-0.04 km s" 1 , HIP 56948 and 
HIP 100963 are essentially solar. 

• AAp e : Regarding the mctallicity, HIP 100963 is al- 
most indiscernible from the Sun 0.01 dex) and 
HIP 56948 (~ +0.01-0.02 dex) is also near-solar (or 
very slightly metal-rich?), while HIP 79672 appears 
to be somewhat metal-rich (~ +0.05 dex). 

• A Li : HIP 79672 and HIP 100963 are markedly over- 
abundant in Li compared to the Sun by ~ +0.7- 
0.8 dex (by a factor of ~ 5-6), while the differ- 
ence from j4li,q is much milder for HIP 56946 (only 
- +0.2 dex or -60%). 

• (v rt ) (equivalent to t> c sini): HIP 56948 has al- 
most the same projected rotational velocity as the 
Sun, while HIP 79672 and HIP 100963 show slightly 
higher values by — 5-10%. 

• A Bc : All three stars (HIP 56948, HIP 79672, and 
HIP 100963) have essentially the same Be abun- 
dances as the Sun, which means that Be does not 
conform to the behavior of Li showing a diversity. 
Consequently, whichever mechanism changing the 
surface Li abundance of these solar twins can not 



influence Be; e.g., if the variation of An is caused 
by an envelope mixing, it should not be so deep as 
to affect Be (see also Randich et al. 2002). 

Taking all these results into consideration, we can draw 
the following conclusions. 

HIP 56948 is surely most similar to the Sun among these 
three stars, not only from the similarity of stellar parame- 
ters but also from the viewpoint of surface Li abundance; 
it may thus deserve the name of "closest ever solar twin." 
However, unlike the argument of Mclendez and Ramirez 
(2007) who derived A,4 Li = -0.02(±0.13), since the at- 
mospheric Li abundance of this star is marginally higher 
than the solar value by ~ 0.2 dex, we still can not call it 
a "genuine" solar twin. Another concern about this star 
is that its luminosity derived from the Hipparcos paral- 
lax appears to somewhat larger than the solar luminosity, 
which in effect makes the age older (cf. table 1). We 
suspect that this inconsistency is attributed to the error 
in the parallax because HIP 56948 is comparatively dis- 
tant. The possibility that its actual n is by — 10% larger 
than the catalogued value may as well be considered, since 
the similarity of T e s, logg, and Ap e should guarantee the 
equality of L (cf. Appendix A of Paper II). 

Regarding HIP 79672 and HIP 100963, they have by 
~ 40 K higher T eS than T Gffi0 , by ~ 0.7-0.8 dex larger A Li 
than j4.Li,0, and by ~ 5-10% larger v e smi than u e sini0. 
Apart from these considerable differences, the parameters 
of HIP 100963 quite resemble the solar values. Meanwhile, 
HIP 56948 shows other noticeable differences from the 
Sun with respect to v t (by +0.03-0.04 km s _1 ) and ^Fe 
(~ +0.05 dex), which makes this star comparatively lower 
ranked as a solar twin among the three. 

4-2. What Controls Lithium? — Roles of Rotation and 

T eff 

Let us turn our attention to the question posed in sec- 
tion 1: "why these solar twins show a diversity in the Li 
line strength in spite of their similarity to one another?" 
Is this attributed to the difference in the rotational ve- 
locity, as suggested in Paper II? According to tables 1 
and 5, the values of v T t(— i; e shn) for the Li-strong group 
(HIP 79672 and HIP 100963) are somewhat larger by ~ 5- 
10% than those for the Li-weak group (HIP 56948 and the 
Sun/Moon), which can also be visually recognized in fig- 
ure 3. Considering that i = 90° for the Sun while i is 
unknown for the three stars, we can assure that the equa- 
torial rotational velocities (v c ) of Li-strong HIP 79672 and 
HIP 100963 are anyhow larger than the solar value (u e ,©), 
which may be just favorable for the working hypothesis of 
Paper II. 

Yet, we feel it still premature to conclude that the stel- 
lar rotation is the only decisive factor to influence the 
surface Li abundance of these solar twin stars. Inspecting 
the core features of Ca n H and K lines of the program 
stars shown in figure 6, which are sensitive to the chro- 
mospheric activity closely related to the stellar rotation 
rate, 6 we see in any of these stars almost no apprecia- 

6 The strength of the Ca II H+K core emission is known to roughly 
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Fig. 6. Emission features in (a) Ca II K and (b) Ca II H 
resonance line cores of the program stars (solid lines), com- 
pared with Kurucz et al.'s (1984) solar flux spectrum (red 
dashed line). The continuum level of each spectrum (with a 
vertical offset of 1.0 relative to the adjacent one) has been 
appropriately adjusted so that the stellar spectrum matches 
the reference solar spectrum at AA ~ ±0.5A. (a) HIP 56948, 
(b) HIP 79672, (c) HIP 100963, and (d) Moon. 

ble differences in comparison to the solar H and K cores 
of Kurucz et al.'s (1984) solar spectrum. (Actually, the 
largest difference relative to the Sun among the four is 
seen in the "Moon" spectrum, which is presumably due 
to the difference in the solar activity phase, because the 
year of 2008 corresponds to the almost minimum activ- 
ity.) Besides, according to Giampapa (2005; cf. figure 
36 therein), the Li-strong HIP 79672 exhibits the stellar 
activity (inferred from Ca II H and K line cores; ampli- 
tude of ~ 10% with period of 8-9 years) quite similar to 
the solar activity. We thus consider that "substantially" 
large difference in the rotational rate is not very likely be- 
tween Li-strong (HIP 79672 and HIP 100963) and Li-weak 
(HIP 56948 and the Sun) groups. 

We rather suspect that the difference in T e g may also 
play a significant role in producing this Li-diversity. As 
reported in Paper II, there is a slanted "lower boundary 
line" at +100 K > AT cff > K in the A Li vs. AT eff 
diagram, below which no stars are found (cf. figure 9 
therein). Interestingly, when we plot HIP 79672 and 
HIP 100963 (both have AT off - +50 K and A Li ~ 1.6) 
on this diagram, they almost fall on this boundary line, 
Au ~ 1 + 1( AT c ff /TOOK), which we regard to be a signif- 
icant fact being worth attention. 

That is, as speculated in Paper II, we consider that the 
diversity of Au (at a given T e g) is due to the difference in 
the rotational velocity (i.e., slower rotators tend to show 



scale with the rotational rate (e.g., Noyes et al. 1984), though its 
exact relationship is still under discussion (see, e..g., Giampapa 
2005). This means, for example, that a solar-type star rotating 
with v c ~ 10 km s _1 would show a stronger core emission by 
several times than the Sun. 



lower Li abundances presumably caused by an enhanced 
envelope-mixing). On the other hand, since the existence 
of steeply-slanted lower boundary of Au means that ev- 
ery slow rotators should settle on this boundary line, the 
Au values of such slow rotators would naturally show the 
marked T cff -dependence of dA u /d(T e g/ 100 K) - 1. This 
scenario may reasonably explain (at least to an order of 
magnitude) the difference in Au by ~ 0.6-0.7 dex between 
HIP 79672/HIP 100963 and Sun/HIP 79672, all show- 
ing superficially slow rotation as the Sun, while T e ff for 
the former being slightly higher than the latter by ~ 40- 
50 K. Consequently, the Li-strong nature of HIP 79672 
and HIP 100963 may be naturally explained by the fact 
that they belong to the "boundary-line stars" (which seem 
to have higher possibilities of hosting planets; cf. subsec- 
tion 5.1 in Paper II). In this sense, we would suggest that 
T e ff is another significant factor (along with rotational ve- 
locity v e ) in controlling the lithium abundances of solar 
twins, especially for slowly-rotating ones. 

Anyway, this is nothing but a phcnomcnological expla- 
nation, and a number of tasks are still left until the real 
physical mechanism involved in determining the surface Li 
abundance is clarified. Hence, investigations (especially 
on the theoretical side) on the inter-relations between ro- 
tation, T e ff, and Au in solar-analog stars are desirably 
awaited, so that the confronted problems could be set- 
tled: 

- Why does such a slanted lower boundary exist in the 
Au vs. T ff diagram of solar-analog stars, below which 
stars do not exist ("forbidden zone")? Any T ff-sensitive 
physical mechanism is acting so as to suppress the further 
Li depiction? 

- Au appears to be positively correlated with both v e sini 
and T e ff. What does this mean? These two factors hap- 
pen to act independently on Au in the same direction? Or 
this is nothing but a superficial effect caused by a tight 
relationship between w c sinz and T e ff? 

- Whichever i> e sinz or T e ff may be the relevant key, the 
physical mechanism working in the envelope of these so- 
lar analog stars must satisfy the condition of changing the 
surface Li without affecting Be. What kind of process is 
that? 

- Finally, from the observational side, the number of 
well-studied solar twins/analogs is still so insufficient as 
to clearly reveal the behavior of Li in Sun-like stars. If 
we could considerably increase the number of the sample 
stars (e.g., ~ 10 3 or even more), it would surely give us 
a new insight to this field (in addition, a nearly-perfect 
solar twin might as well be detected). Besides, given that 
field solar-type stars are diverse in their age (cf. figure 10 
in Paper II), intensively studying the early-G dwarfs in 
old solar-age clusters (e.g., M 67) would also be beneficial 
for disentangling the roles of various stellar parameters on 
this Li problem. 

5. Conclusion 

An intensive spectroscopic study based on the high- 
quality spectra obtained with Subaru/HDS was performed 
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for HIP 56948, HIP 79672, and HIP 100963 (along with 
the Sun/Moon as the reference standard), known to be the 
representative solar twins, in order to (1) clarify which of 
the three most resembles the Sun by precisely establishing 
the various differential parameters relative to the Sun and 
(2) investigate the reason why appreciable differences in 
the surface Li abundance are observed for these superfi- 
cially very similar stars. 

The standard atmospheric parameters (T^ d , logg std , 
vf d , and {Fe/H} std ) were first evaluated by using the 
equivalent widths of Fc I and Fe II lines, the rotational ve- 
locity parameter (v Tt ; which is nearly equivalent to w sini) 
was derived from the line-profile width by eliminating 
the effect of the macroturbulcnce, and the lithium abun- 
dance (^Li) was determined from the Li I doublet line at 
~ 6707.8 A. Further, the differences of atmospheric pa- 
rameters (AToff, Alogg, Avt, and A^4pc) relative to the 
Sun were established by using the method of precision 
differential analysis (Paper I). 

While we could confirm that HIP 79672/100963 have 
appreciably higher Li content by a factor of ~ 5-6 as com- 
pared to Sun/HIP 56948, the Be abundances for all the 
program stars (HIP 56948/79672/100963) turned out to 
be essentially the same as the solar value, which indicates 
that Be is not affected by any mechanism causing the vari- 
ation of Li. 

We found that HIP 56948 is most similar to the Sun 
among the three, not only from the similarity of stellar pa- 
rameters (including rotation) but also from the weakness 
of the Li line (however, Au for this star is still slightly 
larger than Au,@ by ~ 0.2 dex; i.e., not perfectly the 
same). It may thus deserve the name of "closest ever so- 
lar twin." Meanwhile, some remarkable differences from 
the solar parameters are recognized in HIP 79672 and 
HIP 100963, which show somewhat higher T cff (by -40 K) 
considerably larger Au (by — 0.7-0.8 dex) and slightly 
higher rotational velocity (by ~ 5-10%). 

We can see a tendency that the Li-strong HIP 79672 
and HIP 100963 have somewhat larger rotational veloc- 
ity by ~ 5-10% than Li-wcak HIP 56948 and the Sun, 
which is consistent with the suggestion of Paper II that 
Au is closely correlated with the stellar rotational veloc- 
ity. However, it does not seem very likely that a sub- 
stantial difference exists in the rotational velocity between 
these two groups, because no essential differences are seen 
in their chromospheric activities (sensitive to stellar ro- 
tation) inferred from Ca II H+K line cores. We rather 
suspect that the overabundance of Li in HIP 79672 and 
HIP 100963 (by - 0.6-0.7 dex) is attributed to the differ- 
ence in T e g (by — +50 K) relative to the Sun, since these 
two stars fall on the T e ff-sensitive slanted lower boundary 
in the Au vs. T e g distribution as reported in Paper II for 
solar analog stars. However, it is not clear whether and 
how these two factors (rotational velocity and effective 
temperature) are mutually related in affecting the surface 
Li abundance, which remains to be further clarified. 
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Table 1. Target stars and their parameters 



Star 


HIP 


Name 


Sp. type 




, std 

log g 




{Fe/H} atd 


7T 


a 1 tt 


logL 


M 


log age 




£71^6708 




No. 


number 






(K) 


(cm s~ 1 ) 


(km s _1 ) 


(dex) 


(mas) 




(Lq) 


(M ) 


(y) 


(km s- 1 ) 


(mA) 




1 


56948 




G5 


5747.9 


4.409 


0.93 


-0.016 


15.0 


0.05 


+0.13 


0.98 


9.90 


2.18 


3.6 


1.13 


2 


79672 


18 Sco 


Gl V 


5771.7 


4.397 


0.97 


+0.011 


71.3 


0.01 


+0.05 


1.01 


9.70 


2.29 


9.7 


1.60 


3 


100963 




G5 


5760.0 


4.411 


0.93 


-0.040 


35.4 


0.02 


+0.02 


1.00 


9.66 


2.27 


12.1 


1.68 







Sun/Moon 


G2 V 


5737.1 


4.420 


0.95 


-0.036 






(0.00) 


(1.00) 


(9.66) 


2.13 


2.2 


0.91 



Notes: 

In columns 5-8 are listed the Tf^ (effective temperature) , log g btd (logarithmic surface gravity), v^ td (microturbulence) , and {Fe/H} fatd 
(~ Ap c d — 7.50; mctallicity) , which are the "standard" atmospheric parameters spectroscopically determined based on the selected Fe I and Fe II 
lines (see subsubscction 3.1.1 of Paper II). Columns 9 and 10 present the Hipparcos parallaxes (ESA 1997) and their relative errors, while the 
values of luminosity, mass, and age are given in columns 11-13, which were derived from the positions on the theoretical HR diagram with the help 
of stellar evolutionary tracks (see subsection 3.3 of Paper II). The (^ r t) (column 14) is the rotational broadening parameter (nearly equivalent to 
the projected rotational velocity v e s'mi) determined from the widths of a number of Fe I and Fe II lines (cf. subsection 3.2). The equivalent width 
of the Li I 6707.8 doublet and the logarithmic abundance of Li (including the relevant non-LTE correction of +0.07 dex for all four stars) in the 
usual normalization of Ah — 12 arc presented in the last columns 15—16. 
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Table 2. Differential analysis of HIP 56948 relative to the Sun. 





AT 


A log g 


Av t 


AA e T e g 


e v 


e Al 


£A 2 


CA 1 


0"A 2 


Ni 


N 2 


056948 - Sun 


+4.3 


-0.015 


-0.02 


+0.015 5.0 0.010 


0.04 


0.007 


0.008 


0.021 


0.021 


196 


18 


-(Sun -056948) 


+1.3 


-0.027 


-0.01 


+0.012 5.0 0.010 


0.04 


0.006 


0.007 


0.020 


0.018 


189 


16 


(056948 - Sun) 


+2.8 


-0.021 


-0.01 


+0.013 


















AT 


Alogp 


Av t 


[indirect analysis] 

AA UT Ug 


a v 














((056948 -Sun)) 


+17.0 


+0.006 


-0.01 


+0.020 3.8 0.012 


0.01 


0.004 












(via 079672) 


+20.7 


+0.018 


-0.01 


+0.024 

















(via 100963) +13.2 -0.006 +0.00 +0.016 [differences of standard parameters] 

AT = 5747.9 - 5737.1 = +10.8 
Alogg = 4.409 - 4.420 = -0.011 
Av t = 0.93 -0.95 = -0.02 
AA = -0.016 - (-0.036) = +0.020 

Notes. 

The results of differential analyses for the case of (i = 1 and j = 0). The brief description of the data in the table is given 
below while Paper I should be consulted for more detailed explanations. (Note that the effective temperature T e g and the Fe 
abundance Afc are abbreviated as T as A, respectively, in this table 2 along with the following tables 3 and 4.) 

(Direct analysis:) 

— 1st row: the results of (AT_j, Alog<?i_j, Avt-j, AAi-j), the possible errors (e-r, e 9 , and e„) involved in these solutions 
(estimated by the procedure described in subsection 5.2 of Takeda et al. 2002), the root-mean-square errors (eAi, £a 2 ) ° n the 
differential abundances (AAi t i-j and AA2,i-j) from Fe I and Fe II lines corresponding to these ambiguities in atmospheric 
parameters, the standard deviations (cta-i and cta 2 ) around the means of AA\ t i-j and AA^i-j, and the numbers (iVi and N2) 
of the used Fe I and Fe 11 lines. 

— 2nd row: the same as the 1st row, but for the inverse case of j — i; i.e., presented are the parameter differences of ( — AT,_i, 
— Alogijj-i, —Avj-i, and — AAj-i) and the corresponding errors (et, e 9 , e v , £An and £a 2 ). 

— 3rd row: averaged solutions of the parameter differences given in the 1st and 2nd rows; i.e., (ATj), (Aloggij), (Avij), and 
<AA SJ ). 

(Indirect analysis:) 

— The first row gives {{ApiQj)) [equation (15) in Paper I] and (((TAp,i()j)) [equation (16) in Paper I], where p denotes each of 
T, logg, v, and A. 

— Meanwhile, in the following two rows are presented the individual {Api(k)j) values [equation (14) in Paper I] for each 
intermediary star k (from which the ((ApjQj)) and ((o"Ap,i()j)) values in the first row were computed). 

(Inset in the lower-right space:) 

— The (i — j) differences of the standard parameters (T, log g, v, and A) given in table 1 . 



Table 3. Differential analysis of HIP 79672 relative to the Sun.* 





AT 


Alogp 


Av t 


[direct analysis] 
AA e T e g 


e« 




£A 2 


CA 1 


0"A 2 


Ni 


N 2 


079672 - Sun 
-(Sun -079672) 


+48.9 
+48.1 


+0.008 
+0.009 


+0.03 
+0.03 


+0.056 5.0 0.010 
+0.053 5.0 0.010 


0.03 
0.03 


0.006 
0.006 


0.006 
0.006 


0.018 
0.019 


0.016 
0.017 


194 
193 


17 
17 


(079672 - Sun) 


+48.5 


+0.008 


+0.03 


+0.054 


















AT 


A log g 


Av t 


[indirect analysis] 
AA ot o g 


a v 


a a 












((079672 -Sun)) 


+39.3 


-0.013 


+0.04 


+0.048 8.8 0.018 


0.00 


0.004 












(via 056948) 


+30.6 


-0.031 


+0.03 


+0.044 

















(via 100963) +48.1 +0.005 +0.04 +0.052 [differences of standard parameters] 

AT = 5771.7 - 5737.1 = +34.6 
Alogg = 4.397 - 4.420 = -0.023 
A« t = 0.97 -0.95 = +0.02 

AA = 0.011 - (-0.036) = +0.047 



*The results of differential analyses for the case of (i = 2 and j = 0). See the notes in table 2 for the details. 
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Table 4. Differential analysis of HIP 100963 relative to the Sun.* 





AT 


A log g 


Aw t 


\c\ irppf anal vqi til 

AA e T e a 






«A 2 




°A 2 


iVi 


N 2 


100963 - Sun 
-(Sun - 100963) 


+37.8 
+38.6 


+0.021 
+0.026 


+0.01 
-0.01 


+0.003 0.0 0.010 
+0.005 0.0 0.010 


0.04 
0.03 


0.005 
0.005 


0.007 
0.006 


0.019 
0.018 


0.013 
0.014 


194 
188 


18 
17 


(100963 -Sun) 


+38.2 


+0.023 


+0.00 


+0.004 


















AT 


A log 5 


Av t 


[indirect analysis] 

A^4 CJT Ug 


a v 














((100963 -Sun)) 


+33.2 


+0.018 


-0.01 


+0.004 5.4 0.009 


0.00 


0.002 












(via 056948) 


+27.8 


+0.009 


-0.01 


+0.001 

















(via 079672) +38.6 +0.027 -0.01 +0.006 [differences of standard parameters] 

AT = 5760.0 - 5737.1 = +22.9 
Alogfl = 4.411 - 4.420 = -0.009 
Av t = 0.93 -0.95 = -0.02 
AA = -0.040 - (-0.036) = -0.004 

*The results of differential analyses for the case of (i = 3 and j = 0). See the notes in table 2 for the details. 



Table 5. Summary of differential parameters relative to the Sun 



Star 


AT off 


A log g 


Avt 


AAfo 


AA hi 


(«rt)/(«S> 




(K) 


(dex) 


(km s _1 ) 


(dex) 


(dex) 




HIP 56948 


+2.8 


-0.021 


-0.01 


+0.013 


+0.22 


1.02 




+17.0 


+0.006 


-0.01 


+0.020 






HIP 79672 


+48.5 


+0.008 


+0.03 


+0.054 


+0.69 


1.08 




+39.3 


-0.013 


+0.04 


+0.048 






HIP 100963 


+38.2 


+0.023 


+0.00 


+0.004 


+0.77 


1.07 




+33.2 


+0.018 


-0.01 


+0.004 







Notes: 

Columns 2 through 5 present the "star— Sun" differences for each of the parameters {T e s, logp, Vt, and Afo; final averaged results 
extracted from tables 2-4) resulting from the differential analysis based on the method of Paper I, where the values in the upper 
and lower row correspond to the direct solution ((Apio); cf. section 3 in Paper I) and the solution obtained via intermediary 
stars (((ApiQo)); cf. section 4 in Paper I), respectively. Given in columns 6 and 7 are the differential Li abundance relative to 
the Sun and the "star/Sun" ratio of (v r t) (rotational broadening parameter), respectively, which were simply obtained from the 
data in table 1. 



